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1.  INTRODUCTION 

There  are  several  technology  areas  of  interest  to  the  Navy  for  which  an  under- 
standing of  multiple  scattering  of  infrared  radiation  by  fog  is  an  important  considera- 
tion. One  application  is  in  the  field  of  electro-optics  meteorology  (EO-Met),  the 
limitations  of  electro-optical  systems  due  to  water.  Based  on  this  information,  the 
Navy  can  determine  the  utility  of  electro-optical  systems  under  various  field 
conditions . 

An  example  is  optical  communications  between  a satellite  and  a ship.  The  light 
beam  must  traverse  atmospheric  scenes  with  clouds  and  fog.  The  communication 
system  could  use  laser  heterodyne  techniques,  as  discussed  by  Yural^^  Lutomirski 
and  Buserl^^  andJ.P.  Gollub  et  al(^,^  for  example.  However,  since  the  usefulness 
depends  upon  the  laser  beam  coherence,  it  is  essential  to  know  its  degradation  by 
multiple  scattering.  It  would  be  useful  to  correlate  multiple  scattering  data  taken 
at  different  wavelengths  under  different  weather  conditions.  It  is  important  to 
relate  visibility  to  atmospheric  transmittance  at  3.8  and  10.6  laser  radiation 
in  order  to  know  the  propagation  characteristics  of  DF  and  CO2  laser  beams. 

Another  useful  application  is  in  covert  communications  where  beam  blovoip  by 
multiple  scattering  (in  addition  to  atmospheric  turbulence,  such  as  described  by 
Lu'ominski  and  Yura^^^  could  permit  the  opposition  to  intercept  Information  that  is 

Yura,  "Optical  Heterodyne  Signal  Power  Obtained  From  Finite  Sized  Sources 
of  Radiation.  " Appl.  Opt.  13.  150-157  (1974). 

^^^R.F.  Lutomirski  and  R.G.  Buser,  "Mutual  Coherence  Function  of  a Finite  Optical 
Beam  and  Application  to  Coherent  Detection,  " Appl.  Opt.  12.  2153-2160  (1973). 

Gollub  et  al. , "Optical  Heterodyne  Measurement  of  Cloud  Droplet  Size 
Distributions."  Appl.  Opt.  12,  2838-2842  (1973). 

^^^R.F.  Lutomirski  and  H.T.  Yura,  "Propagation  of  a Finite  Optical  Beam  in  an 
Inhomogeneous  Medium,"  Appl.  Opt.  10.  1652-1658  (1971). 
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being  transniitted.  Alternatively,  knowledge  of  multiple  scattering  would  be  useful 
in  applications  where  it  is  desired  to  intercept  or  monitor  the  transmission  of 
Information  by  an  adversary. 

At  present,  the  Navy  resorts  to  theoretical  calculation  for  predicting  multiple 

scattering  of  infrared  radiation  in  fogs  and  clouds,  using  experimentally  determined 

data  on  refractive  index  and  particle  size  distribution.  (See,  for  example,  the 

detailed  IRIA  state-of-the-art  review  performed  for  ONRby  A,  J.  LaTocca  and 

R,E.  Turner;  ')  However,  theoretical  models  of  multiple  scattering  should  be 

checked  with  experimental  results,  especially  in  the  Infrared  beyond  1 >xm  where 

data  is  lacking.  (For  instance,  multiple  scatter  data  on  clouds  and  fog  have  been 

done  at  0.694  nm  radiation  by  S.R.  Pal  and  A. I,  Carswell^®^  and  at  0.9  pm  by 

(7) 

J.V,  Sinstanley  and  C.  Wigmore;  ')  A multiple  scattering  measurement  facility 
w’ould  also  provide  a means  for  comparing  the  accuracy  and  practicality  of  scattering 
theories  and  calculational  methods. 

An  additional  area  of  potential  practical  Importance  is  knowledge  of  the  multiple 
scattering  properties  of  fog  particles  with  different  concentrations  of  salts  from  the 
ocean.  The  salt  could  significantly  affect  beam  scattering,  the  extinction  of  abeam, 
and  the  multiple  SLittering  signal.  (See,  for  example,  single  scattering  data  by 
W.C.  Wells  et  all^^  and  by  H.R.  Calon  et  al(^^  The  method  discussed  in  this  report 
has  the  potential  capability  of  generating  sea- spray-like  fogs. 

^^^A,  J.  LaRocca  and  R.E.  Turner,  ’’Atmospheric  Transmittance  and  Radiance; 
Methods  of  Calculation,”  IRIA  State-of-the-Art  Report,  June  1975,  pp.  30-33, 
54-84,  441-467. 

^®^S.R.  Pal  and  A. I.  Carswell,  ’’Multiple  Scattering  in  Atmospheric  Clouds:  Lidar 
Observations.”  Appl.  Opt.  15.  1990-1995  (1976). 

(7i 

' 'J.V,  Winstanley  and  C.  Wigmore,  ”A  Multiple  Scattering  Correction  For  Lidar 
Systems,”  AGARD  Conf.  Proc.  No.  183  (Oct.  1975),  pp,  27-1  to  27-13. 

^®W.C.  Wells  et  al. , ’’Aerosol  Distributions  in  Maritime  Air  and  Predicted 
Scattering  Coefficients.”  Appl.  Opt.  16.  654-659  (1977), 

(9) 

H.R.  Carlon  et  al. , ’’Infrared  Extinction  Spectra  of  Some  Common  Liquid 
Aerosols,”  Appl.  Opt.  1598-1605  (1977). 
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It  is  well  known  (see  for  instance  articles  by  Zuev  et  all^^  that  the 
Soviets  are  doing  a large  amount  of  research,  both  theoretical  and  experimental,  in 
the  field  of  light  scattering  in  fogs,  hazes,  and  the  atmosphere  (both  single  and 
multiple) . 

In  summary,  because  multiple  scattering  is  difficult  to  predict  theoretically, 
although  it  is  a tractable  problem,  experimental  data  W’ould  be  useful  in  under- 
standing the  phenomenon. 

Section  2 briefly  reviews  field  measurements  of  fog  parameters.  Section  3 
discusses  absorption  and  single  and  multiple  scattering  of  fog.  Section  4 examines 
requirements  on  fog  generation  for  multiple  scattering  applications,  including 
analyses  on  particle  agglomeration  causes.  Section  5 gives  a detailed  discussion 
of  the  Aerodyne  fog  facility,  included  automated  particle  counting  data  acquisition, 
fog- jet  spray  nozzle  operating  parameters,  fog  chamber,  and  unique  particle  counter 
dilution  technique.  Data  analysis  is  provided  in  Section  7,  especially  variation  of 
nozzle  parameters  on  fog  generation  in  relation  to  natural  fog  measurements. 

Section  8 discusses  two  other  main  fog  generation  mechanisms  investigated:  fog 
particle  condensation  on  nuclei  by  either  mixing  with  cold  air  or  by  cooling  with  a 
cold-w’all  heat  exchange.  The  latter  method  appears  to  hold  the  greatest  potential 
for  controllable  generation  of  large  fog  particles  radius) 


^V.E.  Zuev,  "Laser- Light  Transmission  Through  the  Atmosphere,"  Laser 
Monitoring  of  the  Atmosphere,  E.D.  Hinkley,  Ed.  (New  York,  Springer  Verlag, 
1976). 

V.E.  Zuev,  Propagation  of  Visible  and  Infrared  Radiation  In  The  Atmosphere, 
(New  York,  Wiley,  1974). 

'V.E.  Zuev,  M.V.  Kabanov,  andB.A.  Savelov,  "Propagation  of  Laser  Beams 
In  Scattering  Media,"  Appl.  Opt.  8,  137-141  (1971). 
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2.  NATURAL  FOG  MEASUREMENTS  IN  THE  FIELD 


j Table  1 gives  a listing  of  some  measurements  which  have  been  made  on  natural 

fog.  The  average  particle  radius  lies  somewhere  betw'een  0.1  - 10  /im,  and  the  mean 
I or  most  probable  particle  radius  is  somewhere  between  4-6  nm.  The  typical 

f particle  concentration  is  on  the  order  of  100  to  500  cm“3^  with  a maximum  mass 

\ loading  on  the  order  of  0.9  g/m^. 

I Table  1 

I Natural  Fog  Measurements  In  The  Field 


Particle  Radius 
Range  (um) 

Most  Probable 
Particle  Radius 
(^m) 

Maximum 

Concentration 

(em-^ 

Maximum  Mass 
Loading  {g/m3) 

Institute  d'Optlque 
(ArnuK  and  Brlcard)^ 

1 to  8 

4 to  6 

Unknown 

Unknown 

Rand  Corporation 
(Deirmendjian)*^ 

lto9 

4 . 

IOC 

unknown 

Atmospheric  Sciences 
Laboratory  (WSMR) 
(Pinnick,  et  al.)*5 

Ito  10 

5 

Unknown 

0.002  -0.9 

Review  (by  Zuev)  *6 

Unknown 

5 

Unknown 

lO"^  -0.6 

M;  Diem  (Referenced  by 
Yamamoto,  et  al.)*^ 

1 to  3 

4 . 

450 

0.28 

Ai-nulf  and  J.  Erichard,  "Transmission  By  Haze  and  Fog  In  The  Spectral 
Region  0.35  to  10  Microns."  J.  Opt.  Soc.  Am.  47.  491-498  (1957). 


Deirmendjian,  "Scattering  and  Polarization  Properties  of  Water  Clouds  and 
Hazes  In  The  Visible  and  Infrared,"  Appl.  Opt.  3,  187-196  (1964). 

Pinnick,  J.D.  Lindberg,  andE.B.  Stenmark,  "Vertical  Inhomogeneity  In 
Wintertime  Atmospheric  Fog  and  Haze  In  West  Gennany  and  the  Effects  on  IR 
Transmission,"  Unclassified  report  of  the  Atmospheric  Sciences  Laboratory, 
Vliite  Sands  Missile  Range  (1976). 

Zuev,  Propagation  of  Visible  and  Infrared  Radiation  In  The  Atmosphere. 
(New  York,  Wiley,  1974),  pp.  168-177. 

Yamamota,  M.  Tanaka,  and  S.  Asano,  "Radiative  Transfer  In  Water  Clouds 
In  The  Infrared  Region."  J.  Atm.  Sci.  27.  282-292  (1970). 
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A typical  fog  particle  distribution  as  given  by  Deirmendjian^^^^  is: 

N(R)  = 0.024  Njj  R®  e-1-5  R . (2-1) 

Here,  N (R)  is  in  units  of  cin“3  Np  is  in  cm"^,  R is  in  pm.  Figure  la  shows 

a typical  fog  particle  distribution  according  to  Deirmendjian;  note  the  peak  at 
R = 4 pm.  More  realistically,  fog  is  composed  of  a superposition  of  two  or  more 
modes,  as  shown  in  Fig.  lb. 


1 t S 4S<TttlO 

PmrUcIt  ludius  Oio) 


Figure  la.  Typical  Single-Mode  Fog  Particle  Distribution 
(Deirmendjian)! 
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Figure  lb.  Typical  Bimodal  Fog  and  Haze  Particle  Size  Distribution.^ 

The  Circles  Represent  Actual  Measurements.  The  Solid  Line 
Represents  a Bimodal  Curve  Fit  Through  the  Data.  The  Haze 
and  Fog  Modes  Are  Similar  to  "Haze  M"  and  "Cloud"  Models 
of  Delrmendjian. 

Figure  1.  (Continued) 


^ ^C.  Tomasl  and  F.  Tampieri,  "Size  Distribution  Models  of  Small  Water  Droplets 
In  Mist  and  Their  Volume  Extinction  CoefHcients  at  Visible  and  Infrared 
Wavelengths. " Atm.  Environ.  10.  1005-1013  (1976). 


3.  MULTIPLE  SCATTERING 


Multiple  scattering  occurs  when  the  scattering  mean  free  path  A is  much 
less  than  the  optical  path  L,  and  the  scattering  dominates  absorption.  It  is  conven- 
ient to  set  the  optical  depth  of  the  medium,  which  is  the  reciprocal  of  the  scattering 
coefficient  equal  to  , i.e. , 


'^mfp  ^^scat 


(3-1) 


For  example,  for  a monodisperse  distribution  of  aerosols  of  radius  R , 

" ''ext  = ”d  «ecat  ■ (3*2) 

where  Njj  is  equal  to  the  number  of  particles  per  \mit  volume  and  is  the  Mie 

scattering  efficiency. 

In  order  to  understand  the  angular  effect  of  multiple  scattering,  it  is  useful  to 
first  perform  calculations  for  single  scattering.  Figures  2 through  5 show  the 
proportion  of  scattering  due  to  four  different  radii  water  aerosols  as  a function  of 
angle  from  0 to  180°.  The  refractive  index  used  was  m = 1.363  - 0. 003451  for  a 
wavelength  A = 3.8  fim  (interpolated  from  Irvine  and  Pollack^^^V  The  quantity 
(II  + l2)/2  represents  the  average  of  the  proportional  amount  of  light  scattered  for 
vertical  and  horizontal  polarizations,  as  calculated  from  Mie  scattering  theory. 

It  may  be  seen  that  the  back  scattering  at  180°  is  virtually  nonexistant  for  fog 
particles  with  a radius  greater  than  1 fim;  it  is  well  knowm  that  most  single  scattering 
for  larger  particles  is  in  a forward  lobe  near  d = 0°. 


Irvine  and  J.B.  Pollack,  ’’Infrared  Optical  Properties  of  Water  and  Ice 
Spheres.”  Icarus  8.  324-360  (1968). 
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Figure  2.  Single  Scattering  as  a Function  of  Angle  for  Particle 
Radius  R = 1 um 
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Figure  3.  Single  Scattering  as  a Function  of  Angle  for  Particle 
Radius  R = 2 


Figures  6 through  9 are  plots  of  the  relative  angular  scattering  for  what  might 
be  considered  a multiple  scattering  situation.  In  this  case,  the  assumption  has  been 
made  that  as  a photon  scatters  from  the  aerosol,  it  goes  through  a random  walk 
process  subject  to  the  scattering  shown  on  Figs.  2 through  5.  Therefore,  after 
suffering  several  collisions  the  distribution  and  angle  of  the  radiation  take  on  a 
Gaussian  distribution  which  broadens  in  angular  width,  eventually  extending  to  180  . 
In  each  case,  this  has  been  calculated  by  the  following  formula: 


= exp  {-d^/n0^)  . (3“3) 

Here,  0^  is  taken  to  be  the  half  power  point  of  the  quantity  (Ij^  + ^ 

Figs.  6 through  9 represents  the  number  of  scattering  events.  It  may  be  seen,  for 
instance,  that  for  a particle  radius  of  5 fim  (which  is  a most  probable  radius  for  a 
fog  particle)  that  even  after  50  scatters  there  is  almost  no  backscattering  at  180° 
and  that  most  of  the  scattered  light  is  concentrated  between  0 and  90°. 

Further,  even  though  fog  is  almost  nonabsorbing,  after  much  scattering  and 
traversals  through  the  fog  in  this  random  walk  method,  the  laser  beam  absorption 
can  become  appreciable.  Figure  10  shows  the  results  of  absorption  of  DF  laser 
radiation  due  to  multiple  scattering  in  fog.  The  equation  describing  the  curves  in 
Fig.  10  is  as  follows  (for  a given  monodisperse  distribution); 

^ (-I'abs'-*  ' 

From  Eq.  (3-2)  Np  = lAR^Qgcat^mfp^'  therefore,  Eq.  (3-4)  becomes: 


Y ' Qabs/Qscat> 


(3-5) 


Figure  6.  Multiple  Scattering  for  Particle  Radius  R = 1 
n Is  the  Number  of  Scattering  Events 


It  may  be  noticed  that  for  a particle  radius  of  5 ^m  and  n = 50  scatterings,  the  light 


is  attenuated  to  about  10%.  In  fact,  the  value  of  I^/Iq  never  exceeds  30  - 40%  for 
any  particle  radius  for  n = 50  particle  scatterings.  This  information  taken  together 


with  Figs.  6 through  9 suggests  that  virtually  no  li^t  at  DF  laser  radiation  wave- 
lengths will  be  backscattered  at  180°  without  severe  attenuation,  and  that  what  light 


I does  get  backscattered  due  to  multiple  scattering  will  be  due  primarily  to  the  smaller 

t 

‘ radius  particles. 
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Figure  7.  Multiple  Scattering  for  Particle  Radius  R = 2 jim; 
n Is  the  Number  of  Scattering  Events 
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Figure  9.  Multiple  Scattering  for  Particle  Radius  R = 10  (im; 
n Is  the  Number  of  Scattering  Events 
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Figure  10.  Effect  of  Absorption  on  the  Transmission  of  DF  Laser  Radiation  Due  to 
Multiple  Scattering  in  Fog;  Here.the  Loss  in  Transmission  Due  Only  to 
Absorption  After  (n)  Scattering  Events  Is  Included 
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4.  FOG  GENERATION 


Table  2 presents  a summary  of  some  of  the  known  techniques  for  artificially 
generating  fog  for  experiments.  Perhaps  the  most  important  thing  to  note  is  that 
none  of  the  chambers  were  capable  of  generating  a fog  which  could  last  more  than 
20  minutes,  i.e. , the  fog  dissipated  in  time.  It  may  be  seen  at  the  bottom  of  Table  1 
that  use  of  an  aerosol  spray,  which  will  be  continuously  flowing,  has  the  potential 
for  producing  a nondissipating  fog. 

A description  is  given  below  of  the  required  parameters  for  generating  fog  in 
order  to  perform  multiple  scatter  experiments.  The  essential  design  criterion  is 
that  the  density  of  particles  be  sufficient  to  produce  an  optical  depth  which 

is  about  one-tenth  the  dimension  (in  a given  direction)  of  the  fog  chamber. 

4.1  Required  Particle  Density 

Figure  11  shows  a plot  of  number  of  particles  per  cm^  required  to  perform 
multiple  scatter  experiments  for  two  different  mean  free  paths  as  a function  of 
particle  radii,  with  the  assumption  of  monodisperse  particle  distribution.  For 
particle  radii  in  the  4 to  5 >im  radius  region,  it  may  be  seen  that  for  a mean  free  path 
of  about  50  cm,  approximately  10‘*  cm“^  would  be  required  at  X = 3.8  urn;  for  a 
mean  path  of  10  cm,  approximately  5 x 10^  cm"^  would  be  required.  Those  curves 
were  generated  ^rom  Eq.  (3-2). 

Figure  12  shows  the  required  mass  loading  for  the  multiple  scatter  experiment. 
This  was  obtained  from  Fig.  11  by  multiplying  the  particle  density  by  mass  density 
and  particle  volume.  It  may  be  ^jeeh  that  for  particle  radii  on  the  order  of  4 to  5 )im, 
anywhere  from  3 to  30  g/m^  is'required  for  a multiple  scatter  experiment.  It  may 
be  noted  that  the  published  value  mass  loading  from  commercially  available  sonic 
spray  nozzles  is  on  the  order  of  1000  g/m^.  More  will  be  said  later  on  this  value. 


1 


Table  2 

Artificial  Fog  Generation 


L»borntory 

Technique 

Size  of 
Chamber 
(m>) 

Particle 
Radius 
Range  (tim) 

Most 
Probable 
Particle 
Radius  |im) 

Maximum 

Gbocentration 

(Particles/ 

cm^) 

Mass 

Loading 

(g/m3) 

Comments 

Army  Corp*  of 

Engineers  (N.H.) 

(Munis  and  Delaney)''^ 

Cooling  of  steam 
to  -40c 

■ 

1 to  10 

3.5 

260 

Unknown 

Dissipated  In  10 
to  15  minutes 
after  Initial 
Injection 

Army  Corps  of 

Engineers  (N.H.) 

(Munis  and  Oclaney) 

Cooling  of  steam 
to  -«oc 

4 

2 to  12 

6.5 

(lee  fog) 

170 

Unknown 

Dissipated  In  10 
to  15  minutes 

Siberian  Physics 
Technical  Institute  * 

USSR  (Zuev,  et  al.) 

Unknown 

15 

2 to  12 

4to6 

900 

0.6 

Dissipated 

Ae  ronutronic  -Ford 
(Reisman,  et  al.)^^ 

Cooling  of 
saturated  air 

0.4 

Unknown 
(1  to  9 
estimated) 

2to  4 

1 

Unknown 

Unknown 

Dissipates  in  10 
to  20  minutes 

Aerodyne  Laboratory 
Experiment 

Aerosol  Spray 

250 

1 

Unknown 
(1  to  10 
estimated) 

5 

5 X 10* 
(equivalent  1 
mono- 
disperse 
density) 

3-30 

Non  dissipating 

^ Munis  and  A.J.  Delaney,  ’’SimulLaneous  Measurement  of  Laser  Extinction 

in  Warm  Fog  at  Wavelengths  of  0.6328,  1.15  and  10.6  Microns,”  Army  Research 
Report  343,  October  1975,  Army  Corps  of  Engineers,  Hanover,  N J1 . 

Mimis  and  A.J.  Delaney,  "Measurements  of  Laser  Extinction  in  Ice  Fog 
for  Design  of  SEV  Pilotage  System,  ” Army  Research  Report  302,  August  1972, 
Army  Corps  of  Engineers,  Hanover,  N .H. 

Zuev  et  al. , "Attenuation  of  Visible  and  Infrared  Radiation  by  Artificial 
Water  Fogs."  Izv.,  Atm,  and  Ocean.  Phys.  Ser.  509-516  (1965),  trans. 

Reisman,  G.  Gumming,  and  C.  Bartky,  "Comparison  of  Fog  Scattered  Light 
and  Monochromatic  Incoherent  Light,"  Appl.  Opt.  6,  1969-1972  (1967). 
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Mass  Loading  (g/m  ) 


4,2  Agglomeration 

4.2.1  Rate  of  Agglomeration 


If  one  defines  the  mean  time  between  collisions  of  particle  radius  Rg  and 
velocity  V2  in  a cloud  of  particles  of  Rj^,  V^,  and  density  Nj^j  , then 

-1 


■ ^21  ^21  ^DI 


where  = (Ho  + Ri)^  ^ 


21 


(4-1) 


= the  relative  velocity  between 


particles;  For  this  case,  the  rate  of  agglomeration  of  particle  Rj^  in  a cloud  of 
particles  of  density  is  given  by: 


dN 


D1 


N. 


D2 


dt 


(4-2) 


or 


dlogNjji 

dt  " 2 


N. 


D2 


N. 


(4-3) 


D1 


Thus,  the  agglomeration  time  Tj  (for  the  e"^  point  of  agglomeration)  is 


' *^21  ^21  ^D2 


(4-4) 


4.2.2  Turbulence 

Appendix  A gives  a detailed  description  and  analysis  of  the  effects  of  turbulence 
on  agglomeration  of  fog  particles.  The  results  of  that  analysis  show'  that  turbulence 
effects  will  not  be  significant  in  either  the  region  near  the  spray  nozzle  or  far  away 
from  the  nozzle  in  the  fog  chamber.  If  anything,  gravitational  settling  may  possibly 
have  a more  significant  impact. 
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4.2.3  Brownian  Motion 

According  to  the  equipartition  theorem,  the  mean  square  velocity  of  a particle 
of  radius  and  mass  m^  is  given  by: 


<vf>  = 


3kT 

m. 


(4-5) 


where 


mi  = p-|-  TTR? 


(4-6) 


and  pis  the  mass  density,  k is  Boltzmann's  constant,  T is  the  temperature. 

/ V / 2\  1/2 

In  this  case,  using  Eq.  (4-4)  and  the  approximation  / , the 

agglomeration  rate  is  described  by  the  following  equation 


2 


irkT 


1/2 


(R2  +Hi)' 


-3  -31^/2 

«2  - 


(4-7) 


Figure  13  is  a plot  of  Brownian  motion  agglomeration  of  a Deirmendjian  distribution 
for  = 50  cm.  The  agglomeration  time  is  plotted  on  the  verticle  axis,  the 

particle  radius  R.  on  the  horizontal  axis.  It  may  be  seen  that  the  agglomeration  time 

4 

for  a particle  of  R2  = 1 pm  falling  in  a fog  will  be  no  less  than  10  sec.  For  a 

particle  size  of  R„  = 4 pm,  which  is  a maximum  particle  radius  for  many  fogs,  the 
^ 3 

agglomeration  occurs  no  faster  than  10  sec.  With  a transit  time  of  the  fog  through 
the  chamber  of  about  50  sec.  Brownian  motion  will  be  negligible  as  a cause  for 
water  droplet  agglomeration  in  a fog  chamber. 

4.2.4  Gravitational  Settling 

The  equation  for  the  terminal  or  settling  velocity  V is  given  by: 

s 


4 3 

mg  = P-i"  ""R  = 6ffT)RV„ 

6 


(4-8) 
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(sec) 


(Fog  transit  time  through  chamber  • 50  sec) 


R2  ■ 1 Mhi 


R2  ■ 10  urn 


R2  "4  pm 


5 6 7 8 9 10 


Particle  Radius  R^  (pm) 

Figure  13.  Brownian  Motion  Agglomeration  Rate  for  = 50  cm 

where  m is  the  mass  of  the  particle,  g is  acceleration  of  gravity , p is  the  density  of 
water,  17  is  the  viscosity  of  air.  Solving  for  the  settling  velocity  produces: 


V = BP  M.qi 

s 9,j  • 

Empirically  (from  the  Stanford  Research  Institute  Chart  on  the  Characteristics  of 
Particles)  = 6.8  x 10®  R^*®  for  In  cm/sec,  and  R in  cm.  For  this  situation, 
using  Eq.  (4-4),  the  sigglomeration  time  is  described  by  the  following  equation; 


- . 


(4-10) 


1 

! 


I 


( 

I 

1 


6.8x  10 


r; 


1.9 


R 


1.9 


N. 


D2 


Figure  14  is  a plot  of  settling  velocity  vs  particle  radius.  Also  shown  is  a flow 
velocity,  which  is  about  10  times  faster  than  the  particle  settling  velocity  of 
10  cm/sec,  which  is  for  a 5 /zm  radius  particle. 


i 

t!i 

j ; 
* i 

li 


' Particle  Radius  (urn) 

i-  I 

I 

f!  Figure  14.  Gravitational  Settling.  (Source  for  particle  settling  velocity: 

;i  Stanford  Research  Institute  Chart  on  the  Characteristics  of 

^ Particles,  1961.) 
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Figure  15  shows  the  gravitational  agglomeration  time  for  a mean  free  path  of 
50  cm.  The  curves  in  this  figure  describe  the  effect  of  a larger  particle  of  radius 
R2  falling  through  a particle  distribution  described  by  Deirmendjian.  For  instance, 
it  may  be  seen  that  as  the  size  of  particle  R2  approaches  Rj,  the  time  of  agglomeration 
approaches  infinity  since  both  particles  are  falling  at  the  same  rate  and  one  never 
catches  up  to  the  other.  It  turns  out  for  this  situation  that  as  one  goes  to  larger  and 
larger  particle  radii  R2,  one  might  expect  a larger  amoimt  of  agglomeration  due  to 
faster  settling  velocities.  However,  this  is  offset  by  the  fact  that  larger  and  larger 


R2  “4ijm 


R2  • 7Mni 


(Fog  transit  time  through  chamber  - 50  sec) 
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Figure  15.  Gravitational  Agglomeration  Rate  for  = 50  cm 


J 


particle  radii  are  less  frequent  in  number  according  to  a typical  fog  particle 
distribution.  The  point  at  which  a combination  of  the  two  produces  the  greatest 
amount  of  agglomeration  occurs  for  a value  of  Rg  equal  to  about  7 fim.  In  this  case, 
as  shown  in  Fig.  15,  agglomeration  occurs  in  about  400  sec.  Again,  it  is  worth 
pointing  out  the  fog  transit  time  throu^  the  chamber  can  be  made  to  occur  on  the 
order  of  50  sec,  or  almost  10  times  faster  than  this  agglomeration  rate.  Actually, 
the  situation  is  even  better  than  this;  it  says  that  approximately  10%  of  the  fog  has 
agglomerated  only  by  the  time  it  reaches  the  bottom.  At  the  top,  of  course,  there  is 
virtually  no  agglomeration,  and  at  the  middle  there  is  at  most  5%,  so  this  may  be 
an  overly  conservative  estimate  of  the  effective  agglomeration  time  due  to  gravitational 

settling. 


5.  DESCRIPTION  OF  APPARATUS 


prr 


5.1  Concept 

The  key  concept  in  the  design  of  the  fog  chamber  is  to  obtain  a continuous 
flow  of  fog  so  that  steady  state  operation  is  achieved.  The  components  involve  a 
fog  generator  from  which  the  fog  droplets  are  mixed  with  and  carried  away  by  a 
continuous  stream  of  air  to  obtain  a continuous  flow  of  fog.  By  employing  air  of 
relative  humidity  equal  to  1 , no  evaporation  or  condensation  of  the  fog  droplets 
takes  place,  thus  the  original  particle  size  spectrum  of  the  generated  fog  can  be 
preserved. 

The  spray  nozzle  is  used  to  generate  the  fog.  By  adjusting  the  parameters  of 
the  nozzle,  some  control  of  the  size  spectrum  can  be  achieved.  This  vill  be 
described  in  further  detail  in  Subsection  5.2. 

The  overall  schematic  of  the  experiment  is  shown  in  Fig.  16.  The  different 
components  will  be  described  in  the  following  sections. 


Figure  16.  Overall  Schematic  of  Fog  Experiment  : 

] 

I 

I 
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5.2  Sonimist  Spray  Nozzles 


r 


The  fog  is  generated  by  a spray  nozzle  with  the  droplet  size  spectmim  con- 
trolled by  the  characteristic  of  the  nozzle.  Nozzles  from  a line  called  Sonimist, 
manufactured  by  Heat  Systems  Ultrasonics , Inc.  (38  East  Mall,  Plainsview,  NY  11803) 
are  used.  According  to  the  manufacturer,  this  series  of  nozzles  can  give  adjustable 
droplet  sizes  of  1 to  30  lum  (Fig.  17).  The  operating  conditions  of  these  nozzles  are 
shown  in  Fig.  18.  Two  nozzles  of  different  throughput  rates  were  purchased  - the 

Q 

Model  700-1  with  water  flow  rate  of  20  to  100cm  /min  and  the  Model  900-3  with  rate 

3 

of  600  to  900cm  /min.  The  basic  design  of  the  nozzle  is  shown  in  Fig.  19.  The 
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SONIMIST  SPRAY  NOZZLES 

conventional  S^ftAV  nozzles 
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Figure  17.  Operation  of  Sonimist  Spray  Nozzle 

(a)  Schematic  of  nozzle; 

(b)  Particle  size  distribution 

(Source;  Product  information  sheet.  Heat  Systems  Ultrasonics,  Inc.) 
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Figure  18.  Operating  Conditions  of  Various  Sonmist  Spray  Nozzles 
(^urce:  Product  Information  Sheet,  Heat  Systems 
Ultrasonics,  Inc.) 
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Figure  19.  Schematic  of  Sonimist  Nozzle  (Model  700-1) 
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adjustable  parameters  are  air  pressure,  water  flow  rate,  and  gap  setting  (distance 
between  the  air  jet  outlet  and  the  resonator).  The  spray  pattern  is  sensitive  to  the 
gap  setting.  This  is  shown  in  Fig.  20.  According  to  the  manufacturer,  the  droplet 
size  is  a function  of  the  water  flow  rate  and  gap  setting  and  is  rather  independent  of 
the  air  pressure.  The  recommended  air  pressure  for  operation  is  55  to  60  psig. 

Referring  to  Fig.  19  (for  the  700-1  nozzle),  water  enters  through  the  central 
stem  and  goes  out  of  4 small  holes  of  1/32  in.  diameter  at  the  side.  (There  are  ^ 

8 holes  of  the  same  size  for  the  900-3  nozzle.)  Based  on  water  flow  rate  (~0.4cm  /sec) 
the  exit  velocity  of  the  water  jet  is  ~20  cm/sec.  These  streams  of  water  encounter 
a high  velocity  air  jet  coming  out  from  the  orifice  betv'een  the  central  stem  and  the 
outer  jacket.  (The  clearance  is  about  1 mm.)  Since  the  air  pressure  is  ~60  psig, 
the  air  stream  must  be  choked  at  the  orifice,  and  the  air  velocity  coming  out  should 
be  approximately  sonic.  This  jet  of  air  further  interacts  with  the  resonator  to  give 
a violent  oscillating  flow.  The  frequency  is  given  by  the  gap  size  divided  by  the 
sonic  velocity,  and  amounts  to  ~100  kHz.  The  water  stream  interacts  with  this 
oscillating  flow  for  a distance  of  ~1  mm  (dimension  of  the  annular  opening  of  the 


0.112  In. 


Sonimist  (hlodel  700-1) 
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resonator)  and  is  shattered  into  water  droplets.  The  residence  time  t in  this 
shattering  region  is: 

T ~1  mm/20  cm/sec  = 5 ms 


Hence,  the  w’ater  droplets  are  "hit"'  approximately  500  times  before  they  come  out. 
Therefore,  the  Sonimist  is  a very  efficient  liquid  shattering  device. 


5.2.1  The  Fog  Chamber 


The  fog  chamber  is  designed  as  a closed  circuit  wind  tunnel  driven  by  a 
1000  scfm  blower  (Fig.  21).  There  is  also  the  capability  of  throttling  the  blower  to 
give  a low'er  flow  rate.  Because  the  Sonimist  delivers  ~1500  cm^/sec  of  extra  air 
into  the  circuit,  balancing  ports  were  introduced  to  let  this  extra  air  bleed  out.  The 
cross  section  at  the  test  section  is  75  cm  x 75  cm  which  gives  an  unthrottled  flow 
velocity  of  ~80  cm/sec. 

*L-Te-T»l 


Figure  21.  Schematic  of  Fog  Chamber 
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The  Sonimist  is  located  at  the  top  of  the  chamber  and  is  adjusted  so  that  it  gives 
a horizontal  spray  pattern  (Fig.  20).  The  spray  is  mixed  with  the  chamber  air  and 
flows  through  the  test  section.  The  return  path  of  the  air  goes  through  two  filters  to 
remove  the  water  droplets.  By  running  the  air  in  a closed  loop,  it  is  ensured  that 
the  relative  humidity  of  the  air  is  equal  to  1.  Excess  water  is  drained  out  at  the 

bottom  of  the  chamber. 

The  sampling  cone  of  the  PMS  particle  size  spectrometer  is  mounted  on  the 
side  wall  with  provision  of  limited  travel  in  both  the  vertical  and  horizontal  directions 
so  that  different  parts  of  the  fog  chamber  can  be  sampled. 

5.2.2  The  PMS  CSAS-100  Particle  Size  Spectrometer 

The  fog  droplet  size  distribution  is  measured  by  the  particle  size  spectrometer 
model  CSAS-100  (Classical  Scattering  Aerosol  Spectrometer)  manufactured  by 
Particle  Measuring  Systems,  Inc.  (PMS)  (1855,  557th  Court,  Boulder,  Colorado  80301) 
The  apparatus  is  described  in  detail  in  the  literaturel^^^  It  is  briefly  discussed  here 

for  continuity. 

The  spectrometer  measures  forward  scattered  light  at  0.63  urn  wavelength 
collected  from  a volume  defined  by  the  receiving  optics.  The  signal,  received  by  a 
PM  tube  is  sent  to  a pulse  height  analyzer.  The  pulse  height  is  correlated  to  the 
particle  size, and  the  number  of  pulses  corresponds  to  how  many  particles  of  the 
particular  size  are  recorded.  The  saturation  limit  of  the  system  is  reached  when  more 
than  one  particle  is  present  in  the  sampling  volume.  The  CSAS-100  is  limited  to  a 
particle  density  of  ~2  x 10  /cm  . 

The  spectrometer  is  divided  into  4 size  ranges  as  shown  in  Fig.  22.  The  bin 
size  is  constant  within  one  range,  but  is  different  for  different  ranges.  The  ranges 
also  overlap.  Therefore,  the  counts  have  to  be  normalized  by  the  bin  sizes  before 
the  different  ranges  can  be  combined  together. 


G Knolleriberg,  Preprint  Volume,  International  Conference  on  Cloud  Physics, 
July  1976,  Boulder,  Colorado.  PubUshedby  the  Am.  Met.  Soc. , Boston,  Mass. 
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Figure  22.  PMS  - CSAS  Particle  Measuring  System 
(4  ranges  'W'ith  15  bins  each) 

For  low  particle  concentration,  a long  counting  period  is  required  for  good 
statistics.  The  spectrometer  automatic  counting  time  per  range  is  selectable  from 
1 sec  to  15  min. 

5.2.3  Dilution  Scheme 

In  order  to  achieve  multiple  scattering  at  irdrared  wavelengths,  fog  droplet 
desnities  of  10®-10®/cm^  are  necessary.  This  is  a factor  of  10-100  greater  densitj' 
than  the  PMS  spectrometer  can  handle.  Therefore,  the  sampled  fog  should  be 
diluted  with  a known  amount  of  air  so  as  not  to  saturate  the  spectrometer. 


1 


The  requirements  for  dilution  are; 

1.  The  dilution  air  should  not  change  the  size  distribution  of  the  sampled  fog. 
This  requires  the  dilution  air  to  be  precisely  the  same  temperature  as 
the  fog  and  also  be  of  relative  humidity  equal  to  1 so  that  no  evaporation 
or  condensation  take  place. 

2.  The  amount  of  dilution  should  be  precisely  known.  Since  a high  dilution 

is  desired  (~10  to  100),  the  volume  of  fog  sampled  will  be  small  compared 
with  the  volume  of  dilution  air.  This  requirement  indicates  that  a fog/air 
valve  bleed-off  system  is  appropriate. 

The  schematic  of  the  dilution  scheme  is  shown  in  Fig.  23.  A closed  circulation 
loop  driven  by  a compressor  is  employed.  The  amount  of  fog  sampled  is  controlled 

Fog  Chamber  | 


1 


Figure  23.  Schematic  of  the  Particle  Counter  Dilution  Scheme.  Here  the 
Dilution  Factor  Is  Given  by  the  Ratio  Fj^/Fg 
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by  bleeding  an  amount  of  air  Fg  out  of  the  loop.  The  sampled  fog  is  mixed  with  an 
amount  (F^  - Fg)  of  air  which  has  been  equilibrated  in  temperature  with  the  fog 


chamber  through  a heat  exchanger.  At  the  start  of  a run,  Fg  is  turned  to  maximum 
so  as  to  saturate  the  air  in  the  closed  loop  with  moisture.  Any  excess  water  droplets 
are  removed  by  a filter  in  the  loop.  Then  the  desired  Fg  is  set  while  F^  is  adjusted 
to  a fixed  value  of  400  SCFH,  the  flow  rate  at  which  the  PMS  counter  is  calibrated. 

Under  these  settings  the  number  of  particles  the  spectrometer  measures  should  be 
scaled  by  a dilution  ratio: 

Fl 

Actual  count  = measured  count  x -= 

Figure  24  shows  the  performance  of  the  dilution  scheme.  The  fact  that  coimts 
do  not  reach  zero  when  Fg  is  shutoff  is  entirely  due  to  leakage  in  the  plumbing.  This 
leakage  is  estimated  to  be  6 SCFH  from  the  intercept  of  the  figure  and  may  be  used 
as  a correction  offset  to  F2.  i 

I 

The  lower  limit  to  Fg  may  be  discussed  on  three  grounds.  i 

1.  Nonisokinetic  sampling  error:  From  Fig.  25,  for  particle  diameters  less 

I than  20  ^m,  an  acceptable  bound  is  that  the  ratio  of  the  mean  velocity  of 

the  fog  to  the  sampling  velocity  (U^^/Ug)  should  be  less  than  1.5.  At 
present,  the  mean  flow  velocity  at  the  chamber  is  < 1 m/sec  and  the 
sampling  orifice  diameter  is  0.25  in.  Therefore,  F2  should  not  be 
lower  than  6 SCFH  and  F^/Fg  should  be  less  than  66.  However,  one  can 
always  decrease  the  orifice  diameter  to  increase  Ug.  It  is  the  next  two 
considerations  which  set  a lower  limit  to  F2 . 

2.  Leakage:  When  F2  is  small,  the  leakage  in  the  system  may  completely 

overshadow  F2,  and  an  accurate  dilution  ratio  cannot  be  obtained.  I 


3.  Thermal  and  humidity  control  of  the  dilution  air:  When  F2  is  small,  the 
total  amount  of  fog  droplets  is  small.  Therefore,  precise  control  of  the 
dilution  air  temperature  and  humidity  is  very  important,  otherwise  the 
fog  droplets  may  evaporate  or  be  subject  to  further  condensation,  thus 
altering  the  size  distribution . 
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Based  on  the  above  considerations,  it  may  be  seen  from  Fig.  24  that  a dilution 
ratio  of  ~ 20  was  actually  achieved;  by  lowering  F2,  a dilution  ratio  of  ~ 100 

should  be  achievable.  This  would  extend  the  range  of  the  PMS-CSAS  spectrometer 
to  measure  particle  densities  up  to  ~ 2 x 10  Vcm  . 

5,2,4  The  Computer  Interface 

An  RS232C  interface  was  purchased  for  the  PMS  particle  size  spectrometer 
so  that  data  can  be  directly  interfaced  to  the  Prime  400  computer  at  Aerodyne.  Each 
data  record  contains  the  raw  data,  the  settings  on  the  PMS  instrument,  and  the  time 
of  the  day  the  measurement  is  taken.  The  data  are  processed  to  provide  printer 
plots  of  the  size  distribution  at  the  terminal  on  site  and  more  detailed  plots  at  a 
later  time.  Typical  records,  printer  plots,  and  final  plots  of  the  result  are  shown 
in  Figs.  26  through  28. 
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Figure  26.  Typical  Record  From  PMS  CSAS  Particle  Coimter 
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Figure  27.  Printer  Plot  of  Spectrum  From  PMS  CSAS  Particle  Counter  as 
Displayed  on  Computer  Terminal.  The  vertical_axis 
represents  particle  number  density  (^m"*  - cm""^);  the 
horizontal  axis  represents  particle  radius  (urn). 
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Figure  28.  Typical  Plot  of  Fog  Spectrum 


6.  RESULTS  AND  DISCUSSION 


r 


6.1  Data  Collection 

The  parameters  of  the  700-1  Sonimist  were  extensively  varied.  Figure  29  shows 
a map  of  the  operating  gap  sizes  and  water  flow  rates  being  studied.  In  addition,  the 
air  pressure  supply  was  varied  from  30  to  70  psig  at  a water  flow  rate  of 
60  cm^/min  and  0.085  in.  gap  setting.  The  test  results  are  summarized  in 
Figs.  30  through  32  which  are  plots  of  the  particle  size  distributions  obtained  for 
the  different  operating  parameters. 

(Air  Pressure  60  psig) 


u 

c 


be 

c 

ts 

& 

u 


00 

E 


0,140 


0.120 


0.100 


0.080 


0.060 


0 20 


40  60  80  100 


Water  Flow  Rate  (cm^/min) 

Figure  29.  Test  Matrix  of  Gap  Size  versus  Water  Flow  Rate 
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6.2  Test  Results  and  Comparison  with  Fog  Data 

The  test  results  show  that  the  particle  size  spectrum  generated  by  the  Sonimist 
has  a distinctive  bimodal  distribution  which  peaks  at  around  0.8  and  3 fim  radii. 

The  valley  is  at  around  1 nm  radii.  Figure  33  is  an  overlay  of  Figs.  30  through  32. 

By  changing  the  water  flow  rate,  gap  setting,  and  applied  air  pressure,  the  dN/dR 
curves  move  up  and  down  by  more  than  an  order  of  magnitude.  However,  there  is 
little  change  of  the  bimodal  shape  in  the  spectrum. 

The  obtained  spectrum  does  resemble  that  of  natural  fog.  Figure  34  is  the 
result  of  typical  measurement  of  natural  fog  by  Trusty  In  the  0.2  to  0.6  nm  radius 
region,  both  spectrums  have  about  the  same  power  law  dependence. 

In  the  0.6  to  4 ^m  region,  the  Sonimist  spectra  has  peaks  at  3 nm.  Trusty's 
data  also  shows  a slight  peak  at  1 fjm  W'hich  is  followed  by  a flat  region  at  2 nm.  The 
particle  number  density  produced  by  the  Sonimist,  depending  on  the  operating  param- 
eter, is  10^  to  10^  times  denser  than  that  of  Trusty's  natural  fog  data.  Beyond  4 
both  the  Sonimist  and  Trusty's  spectra  take  the  same  form  again  (dN/dR  varies  as 
the  inverse  fifth  to  sixth  power  of  radius). 

Figure  35  shows  the  range  of  optical  depths  produced  by  the  Sonimist  for 
wavelengths  0.63,  1.06,  3.8,  and  10.6  (jim.  The  optical  depth  is  obtained  from 
integration  of  the  product  of  the  spectra  and  the  total  scattering  cross  section.  The 
minimum  optical  depth  of  the  present  system  is  around  3 to  5 m.  This  is 
about  an  order  of  magnitude  larger  than  what  one  would  like  to  achieve  for  a 
laboratory-scale  multiple  scattering  experiment.  A smaller  optical  depth  can  be 
obtained  by  using  a number  of  sonimists  spraying  together.  The  present  system 
is,  however,  suitable  for  studying  the  single  scattering  optical  properties  of  liquid 
droplets. 


Trusty,  Naval  Research  Lab.,  Code  8586,  Washington,  DC  20375. 
Private  communication  (1978). 
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dN/dR--l/CC/MIfRON 


Figure  34.  Natural  Fog  Spectrum  (Ref.  26) 
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From  the  overlay  of  the  various  spectra  obtained  from  the  Sonimist  (Fig.  33) 
we  see  that  although  the  number  density  shifts  up  and  down  with  different  operating 
conditions,  there  is  little  effect  on  the  shape  of  the  spectrum.  Furthermore,  the 
total  volume  of  water  \mder  the  spectrum  [/«®  4/3  IT  dN/dR)j  amounts  only  to  a 

few  percent  of  the  total  amount  of  water  input  to  the  Sonimist.  TTierefore  most  of 


4( 

the  water  must  go  to  some  large  droplets.  It  is  necessary  to  study  the  fundamental 
physics  of  droplet  shattering  in  order  to  know  the  degree  of  controllability  of  the 
particle  size  distribution  in  such  a process. 


The  most  important  parameter  governing  droplet  shattering  is  the  Weber 
number  defined  as: 


. p D 

1 OD 


where 


= density  of  gas  stream 

Uq,  = relative  velocity  between  droplet  and  gas  stream 
<r  = surface  tension  of  the  liquid 
D = diameter  of  the  droplet 


*Since  droplet  volume  ~R^,  only  a small  number  of  large  particles  is  necessary  to 
take  up  the  water.  Typical  values  of  mass  loading  is: 

^ j.  _ water  mass  flow  rate 

mass  oa  ng  - j-g^g  gf  fgg  chamber 

40  gm/min 
5 X 10^  cm^/sec 

= 1.3  X 10“®  gm/cm^ 

The  mass  of  a 100  pm  diameter  droplet  is  0.5  x 10~®  gm.  Therefore,  it  would  only 
require  3 droplets/cm^  to  account  for  most  of  the  water  from  the  Sonimist. 
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The  Weber  number  is  a measure  of  whether  the  aerodynamic  pressure  applied  by 
the  gas  stream  is  able  to  overcome  the  surface  tension  of  the  droplet  and  thus  cause 
shattering.  The  critical  Weber  number  reported  in  the  literature^^^"  varies 
from  3 to  30.  Figure  36  shows  the  shattering  velocity  required  for  the  breakup 
of  different  sizes  water  droplets.  It  is  difficult  to  predict  the  sizes  of  the  small 
droplets  resulting  from  the  breakup  except  that  they  are  much  smaller  than  the  size 
of  the  initial  droplet.  Therefore,  the  spectrum  consists  of  two  parts.  There  are 
the  droplets  which  do  not  breakup  further.  They  peak  at  around  the  diameter  governed 
by  the  critical  Weber  number.  Then  there  are  the  small  droplets  resulting  from  the 
breakup.  They  populate  the  low  end  of  the  spectrum . hi  most  circumstances  it  is 
difficult  to  obtain  a shattering  velocity  higher  than  sonic.  This  corresponds  to  a 
critical  droplet  diameter  of  30  (Fig.  26).  Since  the  Sonimist  operates  at  sonic 
velocity,  the  spectral  peak  of  the  surviving  particles  should  occur  at  this  diameter. 

The  PMS  particle  counter  only  detects  particles  up  to  20  (im  diameter.  Thus  the 
spectrum  observed  in  the  experiment  represents  the  small  droplets  resulting  from 
the  breakup  rather  than  the  surviving  droplets.  Since  these  droplets  are  formed  in 
a random  manner,  the  spectrum  is  insensitive  to  the  adjustment  of  the  Sonimist 
system  parameters , as  we  observed. 


(27)  — 

Y312-131?(19^)^^®’^®  “ Streams  of  Air, '»  Ind.  Eng.  Chem.  43,  6. 

(28) 

"Mechanism  and  Speedy  Breakup  of  Drops,  ” J.  Appl.  Phvs.  30.  11 

(29) 

P.  Simpkins  and  E,  Bales,  "Water  Drop  Response  to  Sudden  Acceleration," 

J_.  Fluid.  Mech.  55,  4,  629-639  (1972). 
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Figure  36.  Shattering  Air  Velocity  Requirement 


6.4  Summary 

The  Sonimist  spray  nozzle  produces  droplets  with  a size  spectrum  that 
resembles  spectra  measurements  of  natural  fog.  The  optical  depth  of  the  present 
artificial  fog  can  be  down  to  3m.  The  spectral  shape  of  the  size  distribution  is 
found  to  be  insensitive  to  the  Sonimist  operating  parameters.  This  is  because  the 
observed  spectrum  is  populated  mainly  by  the  small  droplets  randomly  formed  during 
the  droplet  shattering. 
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7.  OTHER  METHODS  OF  GENERATING  FOG  DROPLETS 
WITH  A CONTROLLABLE  SPECTRUM 

We  studied  other  means  of  generating  fog  droplets  with  a spectrum  controllable 
in  the  0. 1 to  10  ^m  radius  range.  Preliminary  experiments  were  performed  to  alter 
the  size  distribution  of  the  droplets  from  the  Sonimist  by  condensing  steam  onto  them 
and  b>’  evaporating  them  in  hot  air.  Also , grow-th  of  large  particles  from  sodium 
chloride  nuclei  was  carried  out. 

The  generation  of  monodispersed  aerosols  was  reviewed  by  Fuchs  and  Sutagin^^®^ 
and  more  recently  by  Corn  and  Esmenl^^^  Their  review  can  be  summarized  in 
Fig.  37.  Atomization  does  not  provide  particle  radii  smaller  than  20  f^m.  In  the  1 to 
10  ^m  radius  region,  droplets  of  nonvolatile  liquids  had  been  successfully  grown  by  the 
condensation  method.  However,  the  growth  of  water  droplets  is  difficult  because  of 
the  tremendous  latent  heat  of  vaporization  and  high  vapor  pressure. 

Figure  38  shows  the  amount  of  latent  heat  released  from  condensing  various  size 

5 3 

droplets  at  different  densities.  For  example,  to  condense  10  /cm  of  5 fim  radius 
droplets,  the  heat  release  amoimts  to  0. 1 j/cm.  The  heat  capacity  of  air  is 
10~^  j/cm^/°C.  Therefore,  this  heat  release  amounts  to  a ~ 100°C  increase  of  the 
air  temperature.  To  take  out  this  heat,  we  need  to  either  mix  into  the  initial  nuclei- 
carrying  steam /air  mixture  very  cold  air  or  to  p>ass  this  initial  mixture  through  a 
cold-wall  heat  exchanger. 

7.1  Condensation  of  Steam  Onto  Droplets  Produced  by  the  Sonimist 

Preliminary  experiments  were  carried  out  to  grow  the  droplets  produced  by  the 
Sonimist  by  mixing  first  steam  into  the  Sonimist  spray  and  then  cold  air  into  the  spray. 


Fuchs  and  A.  Sutagin,  "Generation  of  Monodisperse  Aerosols.  ” in  Aerosol 
Science.  Ed.  by  C.  N.  Davies,  Academic  Press  (1960).  Ch.  1. 

^^^^M.Com  and  NA  Esmen,  "Aerosol  Generation,"  in  Handbook  of  Aerosols,  ed. 
by  R.  Dennis,  ERDA  (1976),  pp.  9-39. 
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standard  Deviation 


Coodenution  method:  aerosols  of  DOP,  1**,  2‘*,  3*.  4“*,  3"*;  DBP,  6**;  glycerol  7“*;  tetralin,  8”;  ocunoic  add  9'*; 
stearic  acid  10'*’;  linoleoic  acid,  11**;  sodium  Florida,  12**;  13*’;  14'*’;  cadmium,  IS**. 

Atomization  of  liquids:  capillary  tip  in  an  annular  air-jet,  water,  16**,  17”;  rotating  blades,  mineral  oil,  18**;  pulsating 
pressure,  water,  19**;  rotating  disc,  water,  20**,  DBP,  21**;  electrostatic  atomization,  butylbenzoate,  22**. 

Atomization  of  suspensions : latex,  23** ; spores  of  B.  gtobulis,  24'**. 

Dispersion  of  powders : lycopodium,  2S'**. 


Figure  37.  Characteristics  of  Monodisperse  Aerosols  Generated  by 
Different  Methods  (from  Ref.  30) 
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Figure  38.  Latent  Heat  Release  From  Condensing  Various  Size  Water  Droplets 
at  Different  Densities  (Heat  Capacity  of  Air  ~ 10~3  j/cm^/OC) 


An  equilibrium  calculation  on  the  amount  of  water  condensed  is  reported  in  Appendix  B. 
A density  of  10^/cm^  of  5pm  particle  corresponds  to  a mass  loading  of  50  gm/m  . 

The  calculation  (Fig.  B-2)  shows  that  this  can  be  achieved  by  mixing  liquid- nitrogen- 
cooled  gas  with  the  steam/Sonlmist  spray  mixture. 

Figure  39  shows  the  resiilting  droplet  size  spectrum  for  the  following  cases; 

(1)  Sonimist  only  (water  flow  rate  0.5  gm/sec,  air  flow  rate  1. 5 gm/sec); 

(2)  Sonimist  mixed  with  steam  at  0.  2 gm/sec  (Sonimist  parameters  same 
as  (1)) 

(3)  Same  as  (2)  vidth  the  addition  of  1. 2 gm/sec  of  liquid- nitrogen- cooled 
nitrogen  gas. 

The  spectrum  shows  an  order  of  magnitude  increase  of  particles  up  to  the  1 pm 
range  while  the  spectrum  of  the  larger  size  particles  remains  unchanged. 

The  droplet  grov’th  law  in  a diffusion  controlled  process  (cf.  Appendix  B)  may 
be  described  by 

^ ^ . (c  = proportionally  constant) 

Therefore  the  distortion  of  the  particle  size  spectrum  due  to  the  growth  in  size 
of  the  particles  is 

dN(R.  t)  _ rcN(R.t)'| 
dt  " dR  L R J 

Since  N(R)  ~ r”^  for  the  Sonimist  spectrum,  particles  in  the  small  size  range 
of  the  spectrum  will  grow  much  faster  than  those  in  the  large  size  range.  Hence,  we 
observed  the  tilting  towards  a steeper  slope  of  dN/dR  as  more  condensing  water  is 
available.  The  following  Table  computes  the  ratio  ^Q^gerved  Cases  (1) 

and  (2)  where  is  the  spectrum  obtained  from  Case  (1). 


Figure  39.  Spectrum  From  Growth  of  Sonlmist  Droplets  Showing  Effect 
of  Condensation 
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Table  3 

Ratio  of  Observed  to  Calculated  Particle  Growth 


/N  \ AN  . , 

- Al—  I AN  observed 

\2^}  ^observed  - A(N^/2R) 


According  to  the  theory,  should  be  constant.  The  data  is  not 

accurate  enough  (especially  in  regions  of  low  A(N/R))to  determine  the  quantitative 
validity  of  the  theory. 

7.2  Evaporation  of  the  Sonimist  Produced  Droplets 


We  also  tried  mixing  the  output  of  the  Sonimist  with  hot  air  so  as  to  change  the 
spectrum  by  evaporating  the  droplets.  Figure  40  shows  the  resulting  spectrum  from 
mixing  different  amounts  of  hot  air  with  the  Sonimist  spray.  The  process  is  the  re- 
verse of  the  growth  process  with 

dH  c 

s - < 0 (c  = proportional  constant) 

Therefore,  the  small  size  range  of  the  spectrum  shifts  down  faster  than  the  large  size 
range.  At  a hot-air  flow  rate  of  150  SCFH  (compared  with  the  200  SCFH  air  flow  rate 
of  the  Sonimistl,  practically  all  of  the  particles  are  evaporated. 
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7.3  Growth  of  Droplets  by  Heat  Exchanger  Method 


Instead  of  remoxdng  the  latent  heat  of  vaporization  of  the  condensing  w'ater  via 
mixing  w'ith  cold  air,  we  considered  the  option  of  removing  the  heat  by  passing  the 
steam/nuclei/air  mixture  through  a cold  pipe.  In  the  absence  of  condensation  nuclei, 
the  water  will  be  condensing  on  the  wall.  With  the  presence  of  nuclei,  how'ever,  part 
of  the  water  vapor  condenses  on  the  nuclei  which  then  grow  in  size.  The  growth 
process  is  thermally  limited  because  the  only  mechanism  of  heat  removal  is  via  con- 
duction to  the  wall.  This  method  of  growing  nuclei  to  large  droplets  (up  to  27  f^m 
radius)  was  reported  to  have  been  successfully  employed  for  liquids  of  low'  volatility 
(Ref.  30  - 33).  In  Appendix  C w'e  look  into  the  feasibility  of  such  a scheme  for  growing 
w’ater  droplets.  From  a quasi-one-dimensional  calculation  it  seems  feasible  to  grow' 
particles  into  the  5 to  10  ^m  radius  range  by  cooling  the  wall  to  veiy  low  temperature. 

A preliminary  experiment  was  carried  out  to  verify  the  concept.  The  set-up  is 
shown  in  Fig.  41.  Hot  moist  air  is  prepared  by  bubbling  air  through  hot  w'ater  at  about 
25  cm  /sec.  This  air  then  passes  through  a heated  pipe  w'here  the  nuclei  source  is 
located.  The  nuclei  source  (Fig.  42)  consists  of  a grid  of  number  26  nicrome  w'ire, 
heavily  coated  w'ith  sodium  chloride  crystals  and  heated  electrically  up  to  500°C.  The 
nuclei  are  generated  bj'  the  sublimation  and  recondensation  of  the  sodium  chloride. 

This  mixture  of  moist  air  and  nuclei  then  flows  through  a one- inch-diameter  pipe  w'ith 
L/D  equal  to  10.  The  pipe  is  cooled  by  a coolant  jacket.  Ejqaeriments  w'ere  run  both 
W'ith  room-temperature  running  w'ater  and  with  dry- ice- cooled  acetone  as  coolant. 

The  output  particle  spectrum  is  measured  by  the  PMS  particle  counter  using  a vacuum- 
pulled  sampling  line. 

Figures  43  and  44  are  the  spectra  obtained  for  the  tw'o  coolant  temperatures  with 
the  nuclei  source  heated  to  various  temperatures.  Although  the  data  are  too  crude  for 
absolute  quantitative  evaluation, on  a relative  basis  they  do  show  that  there  is  growth  of  a 


(32) 

'J.Burgoyne  and  L. Cohen,  "The  Production  of  Monodisperse  Aerosols  of  Large 
Droplet  Size,"J.  Coll.  Sci.  8 , 364-366  (1953). 

(3^) 

D. Sinclair  and  V . LaMer,  "Light  Scattering  As  a Measure  of  Particle  Size  in 
Aerosols,  "Chem.  Rev.  44.  262-264  (1949). 
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Figure  42.  Sodium  Chloride  Nuclei  Generator 
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Radius  (fim) 

Figure  43.  Droplet  ^ectrum  From  Nuclei  Growth  Using  Room- Temperature 
(T  = 20°C)  Running  Water  as  Heat  Exchanger  Coolant 


1 


significant  amount  of  large  particles.  For  the  room-temperature  water-cooled 
experiment,  the  spectrum  is  populated  up  to  4 to  5 ^m  radii.  For  the  dry-ice-cooled 
case,  there  are  significant  particles  even  up  to  10  jxm  radii.  The  overall  particle 
density  is  approximately  100  times  greater  than  found  with  natural  fog  (Fig.  lb). 


7.4  Summary 

Preliminary  experiments  were  carried  out  to  produce  fog  droplets  of  controllable 
size  spectrum.  They  are  summarized  as  follows: 


Experiment 


Result 


Growth  of  droplets  from  the  Sonimist  by 
condensing  steam  on  to  them  via  cold  air 
mixing. 


Order  of  magnitude  increase  in  the  number 
of  small  particles.  Particles  w’ith  radii 
greater  tiian  2 f^m  remain  unchanged. 


Evaporation  of  the  droplets  from  the 
Sonimist  by  mixing  wi^  hot  dry  air. 


Order  of  magnitude  decrease  in  small 
particles. 


Growth  of  fog  by  condensing  steam  on  to 
NaCl  nuclei  via  a heat  exchanger. 


Significant  amoimt  of  large  particles  (up  to 
5 to  10  fjtm  radii) ; density  is  approximately 
100  times  more  dense  than  in  natural  fog. 


8.  CONCLUSION 


A continuous  flow  fog  chamber  was  built  with  the  fog  generated  by  s.  Sonimist 
spray  nozzle.  The  droplet  size  spectrum  produced  resembles  that  of  natural  fog  and 
is  10  to  100  times  denser.  An  optical  depth  of  3 to  5m  was  obtained  wth  the 
present  set  up,  and  it  is  feasible  to  decrease  the  optical  depth  (i.  e. , increase  the 
fog  particle  density)  by  an  order  of  magnitude  by  using  multiple  spray  nozzles. 

The  spectral  shape  was  found  to  be  insensitive  to  the  operating  parameters  of 
the  Sonimist.  Preliminary  experiments  were  carried  cut  to  produce  fog  with  a var- 
iable size  spectrum.  The  amount  of  small  droplets  produced  by  the  Sonimist  can  be 
increased  or  decreased  by  an  order  of  magnitude  by  mixing  steam  or  hot  air  with  the 
spray. 

Experiments  were  carried  out  to  grow  sodium  chloride  nuclei  in  steam  by  pass- 
ing the  mixture  through  a cooled  pipe.  Preliminary  results  were  very  encouraging  as 
a significant  amount  of  large  particles  (5  to  10  /jm  radii)  were  obtained  by  changing 
the  temperature  of  the  coolant  in  Ihe  heat- exchanger  pipe.  Fog  densities  achieved  with 
these  controllable  , large-particle  spectra  were  up  to  100  times  greater  than  natural 
fog. 
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APPENDIX  A 

TURBULENCE  EFFECTS  ON  DROPLET  COAGULATION 


A simplified  picture  of  the  structure  of  a turbulent  flow  is  as  follows:  large 
scale  velocity  gradients  give  rise  to  flow  instabilities  which  produce  eddies  with  a 
characteristic  scale  comparable  to  that  of  the  gradients,  and  hence,  to  the  dimen- 
sion of  the  duct,  jet,  or  boundary  layer  producing  the  gradients.  The  turbulent  veloc- 
ity u^  characteristic  of  these  layers  is  some  fraction  of  the  mean  flow  velocity;  for 
instance,  in  jets  and  boundary  layers,  the  rms  velocity  fluctuation  is  of  the  order  of 
5 - 20%  of  the  flow  velocity.  Now,  these  large  eddies  are  themselves  unstable  and 
tend  to  break  up  into  smaller  ones  in  times  of  the  order  of  Essentially  in  one 

"turn  around,  " their  energy  being  wholly  transmitted  to  these  smaller  eddies,  since, 
as  long  as  the  Reynolds  number  R^  = puf/p  of  the  eddies  is  large,  viscous  dissipation 
can  be  ignored.  Here,  pis  the  density  of  the  fluid.  There  is  no  inherent  length  scale 
in  the  fluid  (as  long  as  » 1),  and  so  the  pattern  of  edd>'  breakdown  and  energy 
transmission  to  smaller  scales  repeats  itself  continuously.  The  kinetic  energy  per 
unit  voliune  associated  with  eddies  of  size  H and  velocity'  v is  proportional  to  pv  , and 
it  gets  transmitted  dowTi  to  the  next  smaller  size  in  time  on  the  order  of  f/v,  so  the 
rate  of  kinetic  energy  transmission  is: 


(A-1) 


and  this  is  the 
velocity  is: 

V 


same  at  all  scales  below  ( , 

o 


so  that  at  a scale  t the  characteristic 


(A- 2) 


2 2 1/32 

The  viscous  dissipation  ratep(Vv)  then  is  proportional  to  (v/f ) ~ (f  /f  ) = 

-4/3 

jf  ' and  Increases  as  the  eddy  size  decreases.  (Here  p is  the  viscosity.) 

At  some  size  the  eddy  must  be  able  to  absorb  the  energy  E being  transmitted  from 
^o’  ^ given  by: 
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= € = 


3 

f hr 


(A- 3) 


or 


3 
k'^k 


pv.  ’ jP 


= 1 


(A- 4) 


and  since  Vj^  = then 


l/3\^/‘* 

° ^ (Kalmogorov  length) 


(A- 5) 


physically  then,  no  eddies  smaller  than  jPj^  are  to  be  found  in  the  turbulent 
fluid,  since  viscous  effects  are  then  dominant  over  inertial  effects. 

For  the  purposes  of  assessing  the  effect  of  turbulence  on  coagulation  rate,  this 
distance  jPj^  must  be  compared  to  the  "stopping  distance"  or  "relaxation  distance" 
of  a droplet  of  radius  R being  shaken  by  the  velocity  fluctuation.  Consider  the 
drop  equation  of  motion , 


4 3 

-dT  = 67rpR(v-Vjj) 


(A- 6) 


where  is  the  liquid  density  and  Vj^  the  drop  velocity  when  the  fluid  velocity  in  its 
vicinity  is  v.  This  can  be  written  as: 


^rel  dt 


= v - v^ 


(A- 7) 


where  the  relaxation  time  is  : 


rel  9 p 


(A- 8) 
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and  Is  the  time  it  would  take  for  viscosity  to  eliminate  a given  relative  motion  of  the 
drop.  The  characteristic  "shaking  time"  seen  the  drop  is  the  eddy  time  f/v  (after 
eddies  of  size  f)  and  is; 
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v„(f/V 


,2/3, 1/3 

o 


V 


o 


(A- 9) 


Thus,  smaller  eddies  have  a higher  characteristic  "shaking  frequency",  and, 
for  a given  drop  size,  eddies  smaller  than  a certain  size  wll  shake  in  times 
shorter  than  Hence  the  droplets  will  not  be  able  to  follow  the  motion  and  will 

slip  relative  to  the  fluid.  We  find  by  eqiiating  7^,^^  to  : 


or 


^rel 


.2/3  .1/3 
* slip  0 


'’■refo' 


3/2 


slip 


(A- 10) 


(A- 11) 


It  is  clear  that  if  f^^.p  turns  out  to  be  less  than  f j^,  no  eddy  will  exist  which  will  shake 
violently  enough  to  induce  droplet  slip.  The  drops  are  then  locked  to  their  local  fluid 
element,  and  since  no  two  fluids  elements  can  occupy  the  same  position  at  the  same 
time,  turbulence  cannot  influence  the  rate  of  agglomeration.  This  is  true  even  if 
is  less  than  the  mean  distance  between  droplets;  in  this  case  individual  droplets  may 
in  fact  be  caught  in  different  eddies  and  acquire  a corresponding  relative  motion,  but 
still  they  will  "miss"  each  other  as  they  follow  their  individual  fluid  elements  (pro- 
vided the  drop  size  itself  is  not  larger  than  the  eddy  size). 

3 3 

For  water  droplets  in  air  (Pj^  = 1000  kg/m  , p^^^  =1.2  kg/m  , and 
2 X 10^  kg(m-sec)  one  can  calculate  and  ^slip  for  different  u^,  and  R. 

Table  A-1  gives  these  values. 
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Table  A-1 


I 

I 

I 


Calculated  Values  of  f and  -fgjjp  for  Water  Droplets  in  Air 


^o(“) 

u^  (m/sec) 

B(pm) 

i’kOum) 

0. 1 

10 

1 

26 

3.2 

0.1 

10 

2 

26 

26 

0. 1 

10 

10 

26 

3200 

0.1 

1 

1 

147 

0.1 

0.1 

1 

10 

147 

100 

1 

0.  02 

1 

0.  9 X 10"^ 

1 

0.  02 

10 

4900 

0.  09 

The  first  five  lines  in  the  table  refer  to  conditions  typical  in  the  vicinity  of  the 
atomizer  (small  scale,  high-energy  eddies),  while  the  last  two  lines  refer  to  conditions 
in  the  fog  chamber  itself.  It  is  clear  that  the  only  case  for  which  turbulence  would  be 
of  importance  to  agglomeration  is  that  for  the  largest  drops  (10  pm)  at  the  highest 
eddy  velocity  (10  m/sec  and  = 0.  Im).  Certainly,  in  the  chamber  turbulent  effects 
can  be  ignored. 

The  dividing  line  f hs  : 

3/4  \ 3/2 


2Pj^R"v^/9m^ 
r 1/2 


(A- 12) 


which  can  be  rearranged  to: 


A.  = _3_/^ 


xl/2 


VT^^Dj 


^ 0,  074  (R-ater  in  air) 


(A- 13) 
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Therefore,  if  conditions  are  such  as  to  produce  a Kolmogorov  length  much 
longer  than  about  13  drop  radii,  turbulence  agglomeration  can  be  neglected. 

Figure  A-1  summarizes  the  results  of  the  above  analysis. 

The  mathematical  analysis  of  the  spray  nozzle  and  of  agglomeration  in  the  jet 
spray  is  given  below. 

2 

Momentum  Conservation:  27rrhpu  -^exit  ^ 


average  flow’  velocity  yj  u = 


AeAP  1 


47rtanap  r 
turbulent  region  w'idth  h = 2rtana 


4 X 10  ® m 

tano; 

0.1 

AP 

10®  N/m^ 

P 

1,  2 Kg/m® 

’ T* 


h = 0.  2r 


Take  the  scale  of  the  turbulence  (largest  eddies)  as 


f = 0. 4h 
o 


2 

V = 0.  2 \ u (r.  m.  s. 


turbulent  velocity) 


{ = 0. 08r 

V 

0 

0 r 

Then  the  Kolmogorov  length  is  = ( T; 


i l/3\3/4 


u o I 

p y 


for  air 


= 3.4x10-\ 


The  velocity  of  the  Kolmogorov  eddies  is 


V = 

K K 


'"k  “ 


0.049 
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, Velocity  of  Larpost  Eddies,  m/s 


r 


Jet  Radius,  cm 


Figure  A-1.  Turbulent  Jet  Parameters 
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O'  = 0. 1 radian 


Agglomeration  in  the  Jet 

Turbulent  shaking  dominates  over  gravitational  acceleration  if 


Turb,  accel. 


. t K 

K 


if  more  than  9.  8 m/sec  . This  is  true  if 


L2k_L  300  i^m 


and  this  happens  for 


r < 1 m 


For  a given  drop  radius  R,  the  drops  will  slip  relative  to  the  small  (fast) 
eddies  if  their  stokes  relaxation  time 

■^REL.  9 Pliq. 

is  longer  than  the  Kolmogorov  eddy  time 


"^EDDY 
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For  air,  this  implies 


For  the  jet,  drops  10  or  larger  slip  for  r < 0. 4m,  while  drops  1 pim  or 
larger  \\dll  slip  only  for  r < 4 cm. 

Thus,  the  relative  slip  of  1 and  lO/im  drops  occurs  over  a significant  section  of 
the  jet. 

Maximum  loss  rate  of  small  (1  ^m)  drops  to  large  (10  pm)  drops; 


small  _ g“0*  06  ^ 94 

^^small^o 

The  agglomeration  of  water  droplets  in  the  jet  is  small. 
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APPENDIX  B 


CONDENSATION  OF  WATER  VAPOR  BY  MIXING  ^TH  AIR 
OF  LOWER  TEMPERATURE 

Consider  the  mixing  of  a mass,  m^  , of  air  at  temperature,  T^  , with  of 
steam  at  temperature  T2’  The  equilibrium  product  will  be  a mass  of  m^  of  air, 
m^  of  vapor, and  m^  of  condensed  water  all  at  a final  temperature  T.  From  thermal 
balance  we  obtain : 

^ Tl)  ■ ^ ^ * <V-  V <^2-  T)]  + T)’ 


i i 


y 

I 

i 

I 


where 

c^  = heat  capacity  of  air  ; 
c^  = heat  capacity  of  water; 

Cy  = heat  capacity  of  water  vapor  . 

With  the  equilibrium  condition  that  the  vapor  pressure  present  at  T is  given  by  the 
saturated  vapor  pressure,we  can  solve  for  m for  a given  m^/m^  and  Tj^ . The 

result  of  the  calculation  is  sho'vn  in  Fig.  B-1  and  B-2  for  T^  equal  to  20,  0,  -50, 

o 5 3 ^ 

and  -196  C.  A density  of  10  /cm  of  5 pm  particles  amounts  to  a mass  loading  of 

3 

50  gm/m  . Figure  B-2  shows  that  this  is  achievable  bj'  mixing  in  cold  air  at  drj’ 
ice  to  liquid  nitrogen  temperatures. 
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20 


.12 


Figure  B-1.  Plot  of  Mass  Ratios  of  Condensed  Water  Vapor  (Fog  Particles) 
to  Steam  and  to  Air  for  Different  Temperatures  of  Mixing  Air. 
The  Five  Upper  Curves  Give  the  Temperature  of  the  Steam/Air 
Mixture.  The  Three  Lower  Solid  Cuiwes  Are  the  Ratios  for 
Fog/Steam.  The  Five  Lower  Dotted  Curves  are  the  Fog/Air 
Ratios. 


I ' 
i : 
i , 
i ' 
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APPENDIX  C 


PREPABATION  OF  FOG  DROPLETS  BY  COOLING  OF  WATER  VAPOR 
VIA  A HEAT  EXCHANGER:  COMPUTER  SIMULATION 

We  investigated  the  possibility  of  growing  the  small  fog  particles  into  larger 

droplets  (~- 10  /^m  diameter)  by  heat  exchanger  condensation  methods,  i.  e. , the  latent 

heat  of  vaporization  is  removed  via  a heat  exchanger.  This  method  has  successfully 

f33) 

been  used  in  the  past,  noticeably  in  the  LaMer  Sinclair  generator ' ’ and  improved 

versions  of  it,  for  growing  droplets  of  various  low-vapor-pressure  materials.  Gro\ith 
of  svibmicron  size  fog  droplets  were  also  reported.  ' For  the  growth  of  big  droplets, 
however,  the  high  latent  heat  of  vaporization  and  vapor  pressure  of  water  may  present 
a problem.  Consider  the  laminar  flow  of  a vapor  (carried  by  a noncondensible  gas)  in  a 
cold  pipe.  In  the  absence  of  condensation  nuclei,  almost  all  the  vapor  condenses  on 
the  wall.  If  nuclei  are  present,  however,  part  of  the  vapor  condenses  on  them.  These 
nuclei  become  heated  up  by  the  release  of  the  heat  of  vaporization,  and  no  further  con- 
densation takes  place  \mtil  this  heat  is  conducted  away  to  the  wall.  The  process  is 
therefore  thermally  limited.  At  the  same  time,  vapor  is  lost  to  the  wall  by  diffusion 
and  condensation.  The  resultant  droplet  size  is  determined  by  the  competition  betw-een 
the  nuclei  and  the  waU  for  the  vapor. 

The  process  may  be  modelled  by  the  following  equations  (here  the  initial  particle 
size  is  assumed  to  be  large  compared  with  the  mean  free  path): 

Droplet  mass  balance* 


■ar 


— 0!  (n  - nj)v 
r,  v'  oo  d'  m 


(C-1) 


(34) 

G.  Goyen  and  F.  Pidgeon,  "Production  and  Characteristics  of  Submicron 
Water  Aerosols,  "J.  Coll.  Sci.ll,  697-703  (1956). 

* 

Note  that  this  equation  implies  faster  growth  for  small  particles.  Therefore, 
the  particles  tend  to  grow  to  the  same  size. 
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Droplet  thermal  balance 


I ”d  %<=d  dT  * V^d  - ^'■'d  - ^’•d-ar 

Overall  mass  balance 

dn^  _ -N4,rrg  dr^  1 /A  rA  „ "l 

dt  V "Ht”  “v  r Idr  ^ dr  <«/ 

m ' ' 

Overall  thermal  balance 

dT  1 /d  d \ 

fg^g  -af  = »®g<Td  - ^"d  * Kg  - (jj  r T,j 

Vapor  pressure  relation 

v.p.  (T^) 

*^d  ~ iTr^j  • 

where 

n = vapor  molecule  density 

r^  = droplet  radius 

= mass  diffusivity  of  vapor 

V = molecular  volume 

m 

p = density 

C = heat  capacity  per  unit  mass 

K = gas  thermal  conductivity 

g 

k = Boltzman's  constant 

L = latent  heat  of  vaporization  of  vapor 
N = number  density  of  nuclei  = constant 


{C-2) 


(C-3) 


(C-4) 


(C-5) 
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Subscripts 

d = droplet  parameters 

V = vapor  parameters 


oo  = properties  far  from  the  drop 
g = noncondensable  carrier  gas  parameters 

The  condensation  process  may  be  visualized  as  follow’S.  Initially,  say  n = n , 
and  T^=  T^.  Then  as  time  progresses,  T^goes  down  because  of  conduction  to  the 
wall  (last  term  in  Eq.  (C-4)).  This  process  lowers  T^  (2nd  term  in  Eq.  (C-2))  , lowers 
n , (Eq.  C-5). and  starts  the  diffiision  of  n into  the  droplet  (Eq.  C-6)).  This 
diffusion  is  checked  by  the  latent  heat  produced  (last  term  of  Eq.  (C-2)),  and  further 
condensation  is  not  possible  until  this  heat  is  conducted  away.  WTiile  this 
process  is  going  on,  the  wall  acts  as  a sink  for  n (last  term  in  Eq.  (C-3)  , reducing 

OO 

the  amount  of  vapor  available  for  the  droplet  growth. 

Equations  (C-1)  - (C-5)  are  solved  with  the  following  simpli^ing  assumptions. 

It  is  assumed  that  both  n and  (T  - T ) are  parabolic  profiles.  Wliile  the  exact 

00  00  W 

validity  of  the  assumption  is  questionable,  it  does  provide  the  relevant  physics  of 
producing  constant  sink  terms  for  both  n and  T . 

oO  OO 


1 

r dr 


'1 

„2  It  T / 

R \ tt)  w / 


Therefore  a quasi  one  dimensional  treatment  is  possible.  The  effective  pipe  radius  R 
shoxald  be  interpreted  as  a boimdary  layer  thickness  rather  than  the  actual  radius. 

The  following  matrix  of  design  parameters  are  calculated. 

Inlet  condition:  Air  at  50°,  80°C  with  R.H.=1 
Wall  temperature:  293,  223,  77°K 

(corresponds  to  room,  dry  ice,  and 
liquid  nitrogen  temperature) 

Nuclei  density:  10^,  10^,  loVcm^ 

Pipe  Radius:  lin. 
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The  results  of  the  computation  are  shown  in  the  attached  graphs.  Figures  C-1 
to  C-3  are  for  initial  gas  temperature  Tjjj(0)=  50°C,  while  Figures  C-4  to  C-6  are 
for  T (0)=80°C.  Since  the  controlling  time  scale  is  the  mass  and  heat  diffusion  time 
to  the  wall  (t~RVa  ) (assume  Lewis  number  equals  1),  the  pipe  radius  effect  is  to 
approximately  scale  the  time  as  radius  squared.  Therefore,  separate  calculations 
for  different  pipe  radii  are  not  necessary. 

The  plots  of  the  particle  radius  as  a function  of  time  all  show  growth  ~ Vt7 
which  is  consistent  with  the  diffusion  model.  The  growth  time  is  about  1 to  10  sec. 

The  gro\sth  tapers  off  as  the  vapor  is  exhausted.  In  Fig.  C-3  the  droplet  starts  to 
evaporate  after  it  reaches  the  maximum  size.  This  effect  is  solely  due  to  the  assump- 
tion that  there  is  a constant  flux  of  vapor  towards  the  wall  which  acts  as  a sink. 

Figures  C-7  and  C-8  are  contour  plots  of  the  droplet  size  that  can  be  achieved  with 
various  nuclei  density  and  wall  temperatures.  For  an  input  of  80°C  saturated  moist  air 
the  calculation  show's  that  particles  of  5 to  lOj^m  radii  can  be  achieved  over  a wide 
range  of  nuclei  concentration  if  we  cool  the  wall  to  dry  ice  temperature  (223°K). 
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Figure  C-3.  Droplet  Growth  Characteristics  (see  text 
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